Josephson ^-junctions based on structures with complex normal/ferromagnet bilayer 
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We demonstrate that Josephson devices with nontrivial phase difference < (p g < % in the ground state can 
be realized in structures composed from longitudinally oriented normal metal (N) and ferromagnet (F) films 
in the weak link region. Oscillatory coupling across F-layer makes the first harmonic in the current-phase 
relation relatively small, while coupling across N-layer provides negative sign of the second harmonic. To 
derive quantitative criteria for a <p-junction, we have solved two-dimensional boundary-value problem in the 
frame of Usadel equations for overlap and ramp geometries of S-NF-S structures. Our numerical estimates 
show that <p-junctions can be fabricated using up-to-date technology. 

PACS numbers: 74.45,+c, 74.50.+r, 74.78.Fk, 85.25.Cp 
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I. INTRODUCTION 

The relation between supercurrent 1$ across a Josephson 
junction and phase difference (p between the phases of the 
order parameters of superconducting (S) banks is an impor- 
tant characteristic of a Josephson structure 1 ' 2 . In standard 
SIS structures with tunnel type of conductivity of a weak 
link, the current-phase relation (CPR) has the sinusoidal form 
I s (tp) =Asin(p). On the other hand, in SNS or SINIS junc- 
tions with metallic type of conductivity the smaller the tem- 
perature T the larger the deviations from the sin(<p) form 1 and 
Is((p) achieves its maximum at n/2 < (p < %. In SIS junc- 
tions the amplitude B of second harmonic in CPR, 7?sin(2<p), 
is of the second order in transmission coefficient of the tun- 
nel barrier I and therefore is negligibly small for all T. In 
SNS structures the second CPR harmonic is also small in the 
vicinity of critical temperature 7c of superconductors, where 
A ~ (T c — T). At low temperatures T <C 7c, the coefficients A 
and B have comparable magnitudes, thus giving rise to quali- 
tative modifications of CPR shape with decrease of T. 

It is important to note that in all types of junctions discussed 
above the ground state is achieved at <p = 0, since at <p = n a 
junction is at nonequilibrium state. 

The situation changes in Josephson structures involving fer- 
romagnets as weak link materials. The possibility of the so- 
called "ft-state" in SFS Josephson junctions (characterized by 
the negative sign of the critical current 7c) was predicted the- 
oretically and observed experimentally Il2] |29ll . Contrary to 
traditional Josephson structures, in SFS devices it is possible 
to have the ground state (p g = n (so-called 7T -junctions), while 
the <p = corresponds to an unstable situation. It was proven 
experimentally 30,31 that ft -junctions can be used as on-chip %- 
phase shifters or ^-batteries for self -biasing various electronic 
quantum and classical circuits. It was proposed to use self %- 
biasing to decouple quantum circuits from environment or to 
replace conventional inductance and strongly reduce the size 
of an elementary cell 3 - 2 -. 

In some classical and quantum Josephson circuits it is even 
more interesting to create on-chip <p-batteries. They are <p- 
junctions, the structures having phase difference (p g = <p, (0 < 



\<p\ < k) between superconducting electrodes in the ground 
state. The <p-states were first predicted by Mints 33 for the 
case of randomly distributed alternating 0— and it— Joseph- 
son junctions along grain boundaries in high T c cuprates with 
d-wave order parameter symmetry. It was shown later that 
^-junctions can be also realized in the periodic array of 
and % SFS junctions 34 ' 35 . It was demonstrated that depend- 
ing on the length of or % segments in the array, a modulated 
state with the average phase difference (p g can be generated if 
the mismatch length between the segments is small. This (p g 
can take any value within the interval — % < (p g < %. Despite 
strong constraints on parameter spread of individual segments 
estimated in 36 , remarkable progress was recently achieved on 
realization of ^-junctions in such arrays 37 . 

In general, in order to implement a ^-junction one has use 
a Josephson junction having non-sinusoidal current-phase re- 
lation, which, at least, can be described by a sum of two terms 



7 s (<p) =Asin(q>) +Bsin(2q>). 



(1) 



Moreover, the following special relationship between the am- 
plitudes of the CPR harmonics, A, and, B, is needed for exis- 
tence of equilibrium stable state 38 ^ 



|7?|>|A|/2,7?<0. 



(2) 



In conventional junctions, the magnitude of A is larger than 
that of B and the inequalities (O are difficult to fulfill. How- 
ever, in SFS junctions in the vicinity of to % transition the 
amplitude of first harmonic in CPR is close to zero, thus 
opening an opportunity for making a <p— battery, if B can 
be made negative. It is well-known that SFS junctions with 
metallic type of conductivity, as well as SIFS structure s 40 ' 41 
with high transparencies of SF interfaces have complex decay 
length of superconducting correlations induced into F-layer 
%H = %\ + Unfortunately, the conditions (O are violated 
in these types junctions since the A ~ exp{ — L /|j } cos {L/tp) , 
B~ -exp{-2L/^i}cos(2L/^ 2 ), and forL = {n/2)t, 2 corre- 
sponding to the first 0-7T transition the second harmonic am- 
plitude B is positive. 

Quantitative calculations made in the framework of micro- 
scopic theor y 42 ' 43 confirm the above qualitative analysis. In 
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FIG. 1 : a) The S-NF- S junction, b) the SN - FN - NS junction. 



Ref i ' it was demonstrated that in SFS sandwiches with ei- 
ther clean or dirty ferromagnetic metal interlayer the transi- 
tion from to % state is of the first order, that is B > at any 
transition point 3 -. 

It was suggested recently in— ~— to fabricate the "current 
in plane" SFS devices having the weak link region consisting 
from NF or FNF multilayers with the supercurrent flowing 
parallel to FN interfaces. In these structures, superconductiv- 
ity is induced from the S banks into the normal (N) film, while 
F films serves as a source of spin polarized electrons, which 
diffuse from F to N layer thus providing an effective exchange 
field in a weak link. Its strength it can be controlled 5 2^i by 
transparencies of NF interfaces, as well as by the products of 
densities of states at the Fermi level, Nf, Nn, and film thick- 
nesses, dF, dN ■ It was shown in^dr— that the reduction of 
effective exchange energy in a weak link permits to increase 
the decay length from the scale of the order of ~ 1 nm up to 
^100 nm. The calculations performed in these papers did 
not go beyond linear approximation in which the amplitude 
of the second harmonic in the CPR is small. Therefore, the 
question of the feasibility of <p— contacts in these structures 
has not been studied and remains open to date. 

The purpose of this paper is to demonstrate that the same 
"current in plane" devices (see Fig. [TJ) can be used as effec- 
tive (p-shifters. The structure of the paper is the following. In 
Sec|n]we present general qualitative discussion of the micro- 
scopic mechanisms leading to formation of higher harmonics 
in the CPR. In SecJTTT] we formulate quantitative approach in 
terms of Usadel equations. In Sec [IV] the criteria of <p-state 
existence are derived for ramp-type S-FN-S structure. Section 
[V] shows the advantage of the other geometries in order to re- 
alize <p-state. Finally in SeclVIIwe consider properties of real 
materials and estimate the possibility to realize ^-states using 
up-to-date technology. 



II. CPR FORMATION MECHANISMS 

In this section we shall discuss microscopic processes 
which contribute to formation of CPR in Josephson junctions. 
The physical reason leading to the sign reversal of the coeffi- 
cient B in SFS junctions compared to that in SNS structures 
can be understood from simple diagram shown in Figj2] il- 
lustrating the mechanisms of supercurrent transfer in double 



barrier Josephson junctions. 

Consider electron-like quasiparticle e~ propagating across 
SINIS structure towards the right electrode. This quasiparticle 
can be reflected either in the Andreev or in the normal channel. 

The result of the first process (see Figf2^) is generation 
in the weak link region (with an amplitude proportional to 
exp(ixi)) of the hole h + propagating in the opposite direc- 
tion. Andreev reflection of this hole at the second inter- 
face (with an amplitude proportional to exp(— iXi)) results 
in transfer of a Cooper pair from the left to the right elec- 
trode with the rate proportional to the net coefficient of An- 
dreev reflection processes 52 53 at both SN interfaces, AR((p) — 
a(<p)exp(/<p),<p = (Xz — Xi)- The amplitude, a(<p), depends 
on geometry of a structure and on material parameters. Note 
that for given values of these parameters a((p) = a(—(p), ac- 
cording to the detailed balance relations 52 . Similar consider- 
ations show that a quasiparticle e~ moving towards the left 
electrode generates a Cooper pair propagating from the right 
to the left interface with the rate proportional to AR(—(p) = 
a(<p)exp(— /<p). The difference between two processes de- 
scribed above determines a supercurrent 1$, which is propor- 
tional to sin (<p). 

The result of the second process is the change (with an am- 
plitude proportional to exp(ixi)) of the <?~ propagation direc- 
tion to the left electrode and nucleation of a Cooper pair and a 
hole propagating to the right electrode (with an amplitude pro- 
portional to exp(— i(p)). After normal reflection from the right 
interface (with an amplitude proportional to explixi)) the hole 
arrives at the left SN interface and closes this Andreev loop 
by generating a Cooper pair in the left electrode and an elec- 
tronic state (with an amplitude proportional to exp(— iX\))- 
The Cooper pair have to undergo a full reflection at SN in- 
terface, thus again a pair is generated moving in the direction 
opposite to that in the main Andreev loop. The net coefficient 
of this Andreev reflection process is BR(q>) — j3(<p)exp(2/<p). 
For a quasiparticle e~ moving in the weak link towards the 
left electrode the same consideration leads to generation of 
two Cooper pairs moving from the left to the right with the 
rate proportional to BR(—q>) = j3(<p)exp(— 2/<p). The differ- 
ence between these two processes determines a part of super- 
current Ig proportional to sin(2<p). 

We have shown that supercurrent components proportional 
to sin(<p) and sin(2<p) have opposite signs, and the coefficient 
B in Eq.([T]) is negative. This statement is in a full agreement 
with calculations of the CPR performed in the frame of mi- 
croscopic theory of superconductivityii 2 -. It is valid if a super- 
current across a junction does not suppress superconductivity 
in S electrodes in the vicinity of SN interfaces^"—. In addi- 
tion, an effective path of the particles in the second process 
discussed above is two times larger than in the first one. This 
leads to stronger decay of the second harmonic amplitude B 
with increasing the distance L. 

In SFS junctions the situation becomes more complicated. 
The exchange field, H, in the weak link removes the spin de- 
generacy of quasiparticles. As a result, one has to consider 
four types of Andreev's loops instead of two loops discussed 
above. One should also take into account the fact that wave 
function of a quasiparticle propagating through the weak link 
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FIG. 2: Diagrams of the processes forming the first (a) and second (b) 
harmonics of the CPR in the SNS and SFS structures. 



acquires an additional phase shift <p// proportional to the mag- 
nitude of the exchange field 57 . The sign of <p// depends on mu- 
tual orientations between magnetization of the ferromagnetic 
film and the spin of a quasiparticle. Taking into account these 
phase shifts and repeating arguments similar to given above, 
one can show that the coefficients A and B in Eq.(Q]) acquire 
additional factors cos(2<p#) and cos(4<p#), respectively. At 
the point of "0" - "%" transition the coefficient A = 0, that is 
q>H = 7C/4. As a result, cos (4<p//) provides an additional fac- 
tor, which changes the sign of the second harmonic amplitude 
B in SFS structures from negative to positive. 

In the present study we will show that contrary to SFS 
devices with standard geometry, it's possible to realize q>- 
junctions in the structures shown in Fig. Q] Qualitatively, 
these structures are superpositions of parallel SNS and SFS- 
channels, where supercurrent Is((p) can be decomposed into 
two parts, /#(<?>) and If(<P), flowing across N and F films, 
respectively. For L -C <jj# and at sufficiently low tempera- 
tures /jv(<p) has large negative second CPR harmonic B^. For 
L > |i supercurrent in the SFS-channel exhibits damped oscil- 
lations as a function of L. In this regime the second harmonic 
of CPR is negligibly small compared to the first one. Large 
difference between decay lengths of superconducting correla- 
tions in N and F-materials allows one to enter the regime when 
<jji < L < In this case the first CPR harmonic A = An +Ap 
can be made small enough due to negative sign of Af, while 
the second CPR harmonic B ps Bn is negative, thus making it 
possible to fulfill the condition ©. Note that we are consid- 
ering here the regime of finite interface transparencies, when 
higher order harmonics decay fast with the harmonic order. 
Therefore, it is sufficient to consider only the first and the sec- 
ond harmonics of the CPR in all our subsequent discussions. 

We show below that the mechanism described above indeed 
works in the considered S-FN-S junctions, and we estimate 
corresponding parameter range when <p— states can be real- 
ized. 



III. MODEL 

We consider two types of symmetric multilayered struc- 
tures shown schematically on Fig[T] The structures consist 
of a superconducting (S) electrode contacting either the end- 
wall of a FN bilayer (ramp type junctions) or the surface of 
F or N films (overlap junction geometry). The FN bilayer 
consists of ferromagnetic (F) film and normal metal (N) hav- 
ing a thickness dp, and respectively. We suppose that the 
conditions of a dirty limit are fulfilled for all metals and that 
effective electron-phonon coupling constant is zero in F and 
N films. For simplicity we assume that the parameters Ybn 
and Ybf which characterize the transparencies of NS and FS 
interfaces are large enough 



Ybn - 
Ybf - 



Pa"5jV Pn%n ' 
R BF A BF „ psZs 
Pf^f ^ Pf^f ' 



(3) 



in order to neglect suppression of superconductivity in S parts 
of the junctions. Here Rbn,Rbf and Abn,-Abf are the resis- 
tances and areas of the SN and SF interfaces, §v and %p 
are the decay lengths of S, N, F materials and p$, p^ and pp 
are their resistivities. 

Under the above conditions the problem of calculation of 
the supercurrent in the structures reduces to solution of the set 
of Usadel equation s 3 ' 4 ' 58 



£-5 [G 2 m d^ m ] - ^-<t> m =0,G m = 



ft) 



v / (Q 2 + <t>< <i>* 



(4) 



where <i> m and G ffl are Usadel Green's functions in <t> 
parametrization. They are <P a ^ and Go>,n or <Pa, j? and G^.f 
in N and F films correspondingly, ft) = nT(2m + 1) are Mat- 
subara frequencies (m=0,l,2,...), ft) = ft) + iH, H, is exchange 
field of ferromagnetic material, <i; 2 = ^ N F — Dn f /2nTc for 
N and F layers respectively, Dn.f are diffusion coefficients, 
d = (d /dx,d /dz) is 2D gradient operator. To write equations 
©, we have chosen the z and x axis in the directions, respec- 
tively, perpendicular and parallel to the plane of N film and 
we have set the origin in the middle of structure at the free 
interface of F-film (see FigQ]). 

The supercurrent Is((p) can be calculated by integrating the 
standard expressions for the current density jN,F(<P,z) over 
the junction cross-section: 



2ejN.F (<PJ 

TlT 



dp 

W J j F {(p,z)dz^ 





W f j N ((p,z)Wdz, 
d F 



(5) 



where W is the width of the junctions, which is supposed to 
be small compared to Josephson penetration depth. It is con- 
venient to perform the integration in (0 in F and N layers 
separately along the line located at x = 0, where z-component 
of supercurrent density vanishes by symmetry. 

Eq.© must be supplemented by the boundary conditions 59 . 
Since these conditions link the Usadel Green's functions cor- 
responding to the same Matsubara frequency ft), we may sim- 
plify the notations by omitting the subscript ft). At the NF 
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interface the boundary conditions have the form: 



HK*&r^ = g(**-S**) 



R-BFN-A-BFN PF^F 

Ybfn — — z — E = Ybnf - 



(6) 



where Rbfn and Abfn are the resistance and area of the NF 
interface. 

The conditions at free interfaces are 



d<P N = Q d<Z> F = Q 



dn 



dn 



(7) 



The partial derivatives in (0 are taken in the direction normal 
to the boundary, so that n can be either z or x depending on the 
particular geometry of the structure. 

In writing the boundary conditions at the interface with a 
superconductor, we must take into account the fact that in our 
model we have ignored the suppression of superconductivity 
in electrodes, so that in superconductor 



* s (±L/2)=Aexp(±ip/2), G s 



CO 



\Zco 2 + A 2 ' 



(8) 



where A is magnitude of the order parameter in S banks. 
Therefore for NS and FS interfaces we may write: 



YbnZn^- = ^ <&jv - *s(±L/2)) , 
an Gn 



(9a) 



Ybf^f^- = ^-Uf- -* s (±L/2)) . (9b) 
dn Gf \ CO J 

As in Eq. (O, n in Eqs. d9al i, j9bl is a normal vector directed 
into material marked at derivative. 

For the structure presented in FigQJ, the boundary- value 
problem (01 - d9bl was solved analytically in the linear 
approximation 47,48 , i.e. under conditions 



G N = sgn(co), G F = sgn(co). 



(10) 



In the present study we will go beyond linear approximation 
where qualitatively new effects are found. 



IV. RAMP-TYPE GEOMETRY 

The ramp type Josephson junction has simplest geometry 
among the structures shown in Fig[T] It consists of the NF 
bilayer, laterally connected with superconducting electrodes 
(see FigOJ). 

In general, there are three characteristic decay lengths in the 
considered structure 44 - ' 47 ' 68 . They are §v, = §j +1^2, and 
£ = £1 + iQi- The first two lengths determine the decay and 
oscillations of superconducting correlations far from FN in- 
terface, while the last one describes their behavior in its vicin- 
ity. Similar length scale £ occurs in a vicinity of a domain 



wall 6 -^— . In the latter, exchange field is averaged out for an- 
tiparallel directions of magnetizations, and the decay length of 
superconducting correlations becomes close to At FN in- 
terface, the flow of spin-polarized electrons from F to N metal 
and reverse flow of unpolarized electrons from N to F sup- 
presses the exchange field in its vicinity to a value smaller than 
that in a bulk ferromagnetic material thus providing the exis- 
tence of £. Under certain set of parameters 44 these lengths, 
£1, and, £2, can become comparable to which is typically 
much larger than §1 and §2, which are equal to E,f y/ TtTc / H 
for H > kT c . 

The existence of three decay lengths, £, and <ij#, should 
lead to appearance of three contributions to total supercurrent, 
^Vj I fn and Ip, respectively. The main contribution to 1^ com- 
ponent comes from a part of the supercurrent uniformly dis- 
tributed in a normal film. In accordance with the qualitative 
analysis carried out in Section II, it is the only current compo- 
nent which provides a negative value of the amplitude of the 
second harmonic B in the current-phase relation. The smaller 
the distance between electrodes L, the larger this contribution. 
To realize a <p— contact, one must compensate for the ampli- 
tude of the first harmonic, A, in a total current to a value that 
satisfies the requirement (fj). Contribution to A from also 
increases with decreasing L. Obviously, it's difficult to sup- 
press the coefficient A due to the Ifn contribution only, since 
Ifn flows through thin near-boundary layer. Therefore, strong 
reduction of A required to satisfy the inequality (f2]i can only 
be achieved as a result of compensation of the currents In and 
If flowing in opposite directions in N and F films far from FN 
interface. Note that the oscillatory nature of the Ip (L) depen- 
dence allows to satisfy requirement (0 in a certain range of L. 
The role of Ifn in a balance between and Ip can be under- 
stood by solving the boundary value problem (0]) - d9bb which 
admits an analytic solution in some limiting cases. 



A. Limit of small L. 

Solution of the boundary-value problem (|4"l)-(|9bb can be 
simplified in the limit of small distance between supercon- 
ducting electrodes 



L<min{4i,^ v }. 



(11) 



In this case one can neglect non-gradient terms in dU and ob- 
tain that contributions to the total current resulting from the 
redistribution of currents near the FN interface cancel each 
other leading to Ifn = (see Appendix lAl for the details). As 
a result, the total current Is(cp) is a sum of two terms only 

I S ((P)=In{(P)+If{(P), 



1 



2eI N (<p) = 

KTWd N Ybn^nPn 



A 2 GnGs sin(<p) 



co- 



1 



2eI F (q>) = 

KTWdp Ybf&Pf 



A 2 GnG s sin(<p) 



co- 



in) 



(13) 
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where Gn — , , % , a ■ The currents /*/(<?) and Zf(©) 

flow independently across F and N parts of the weak link. The 
Inf(<P) dependencies coincide with those calculated previ- 
ously for double-barrier junctions 59 in the case when L lies 
within the interval defined by the inequalities ( fTTT i. 

It follows from dT2l . ( TT31 that in the considered limit nei- 
ther the presence of a sharp FN boundary in the weak link 
region, nor strong difference in transparencies of SN and SF 
interfaces lead to intermixing of the supercurrents flowing in 
the F and N channels. It is also seen that amplitude of the first 
harmonic of If{<p) current component is always positive and 
the requirement (f2]l can not be achieved. 

B. Limit of intermediate L. 



C. <p-state existence 

The conditions for the implementation of a (j?— contact are 
the better, the larger the relative amplitude of the second har- 
monic which increases at low temperatures. Therefore, low 
temperature regime is most favorable for a <p— state. In the 
limit T <C Tc we can go from summation to integration over 
CO in (Q2K (H3, ( iBBT l- dBTTI i. From (HU we have 



2eI N ((p) A , e> N 

-W7^ = v P n * sin% sin 9 , 
Wd N Ybn^nPn 2 



(19) 



where K(x) is the complete elliptic integral of the first kind. 
Expanding expression dT9b in the Fourier series it is easy to 
obtain 



For intermediate values of spacing between the S electrodes 



(14) 



and for the values of suppression parameters at SN and SF in- 
terfaces satisfying the conditions (01, the boundary problem 
(|4l-(|9bl can be solved analytically for sufficiently large mag- 
nitude of suppression parameter Jbfn- It is shown in Appendix 
IE1 that under these restrictions in the first approximation we 
can neglect the suppression of superconductivity in the N film 
due to proximity with the F layer and find that 

<% = Acos(-) + * 2 j- 

2 Ybn<jn Qv 



Gn — 



CO 



co 2 + A 2 cos 2 (^) 
(15) 

while spatial distribution of <S>f(x,z) includes three terms: the 
first two describe the influence of the N film, while the last 
one has the form well known for SFS junctions^— ^. 

Substitution of these solutions into expression for the su- 
percurrent (01 leads to Is((p) dependence consisting of three 
terms 



Is(<P) = M<P) +M<P) +Ifn(<P)- 



(16) 



Here /jv(<p) is the supercurrent across the N layer. In the con- 
sidered approximation /v(<p) is given by the expression (fT2l) . 
The second term in ( TToT l equals to supercurrent across SFS 
double barrier structure in the limit of small transparencies of 
SF interfaces 69 .^ 



2elp((p) A 2 sin(<p) 
KTWd F jIf^fPf 



• CO 2 V H sinh (2qp) 



(17) 



where q L = Ly/o./2£ F , Q. = |Q| + iH sgn(Q.)/7tT c , Q = 
co/kTc. 

The last contribution is shown in|B]to contain three compo- 
nents 



Ifn(<P) 



*FN1 



((P)+Ifn2{(P)+Ifn3((P)- 



(18) 



with additional smallness parameters Ybfn an( ^ Ybfn^ f /^ n 
compared to the current If{<?) given by Eq.dT7T>. Neverthe- 
less, these currents should be taken into account in the anal- 
ysis because they decay significantly slower than Ip{(p) with 
increasing L. 



A N = Qo- 
n 



B N = 2A N Qq I xWl-x 2 K(x)dx = T B Q , (21) 

n J 
o 

where Qq = AWdN/ey BN ^ N p N , A N , B N are the first and the 
second harmonic amplitudes of In((P), 



Wl-x 2 K(x)dx = Y A Q , 



(20) 



T A = 



2n z 



r 2 (-^)r 2 (|) 



0.973, 



T B = 2T A -| 3 F 2 Q,i |;1,4;1 ] - -0.146. 

where F(z) is Gamma-function and p F q is generalized hyper- 
geometric function. 

Evaluation of the sums in ( fTTT i, (IB 1 51 >- (IB 17b can be done 
for// » nT c and T < T c resulting in I F {<p) =A F sin (<p) with 

A F =P -^=exp(-KL) cos (kL + ^ , (22) 

K = Vh/ V2&, h = H /nT c andPo = AWd F /eYBF&PF- Sub- 
stitution of (l20t . (LTTh into the inequalities (01 gives <p-state 
requirements for ramp-type structure 



e 



^9 IT I ~ _ ^VBF d N ^FpF 

<2 I B , £- — (23) 

2jBN d F C,NPN 



*F(L) = exp(-KL)cos( fcL+- 



This expression gives the limitation on geometrical and ma- 
terials parameters of the considered structures providing the 
existence of ^-junction. Function ^(L) has the first mini- 
mum at kL — 7t/2, *P(7t/2jc) m —0.147. For large values of 
e inequality (l23l can not be fulfilled at any length L. Thus 
solutions exist only in the area with upper limit 



e <= y( *ffi« 0.216. 



(24) 
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At e « 0.216 the left hand side of inequality (1231 ) equals to its 
right hand part providing the nucleation of an interval of kL 
in which we can expect the formation of a ^-contact. This in- 
terval increases with decrease of £ and achieves its maximum 
length 



1.00 < kL< 2.52, 



(25) 



at e = y a +2 \ t b \ ~ 0-116. It is necessary to note that at £ = 
— 1 P(^/2k - )/Ta 0.151 there is a transformation of the left 
hand side local minimum in (23) , which occurs at kL = n/2, 
into local maximum; so that at £ rj 0.116 the both sides of 
d23l become equal to each other, and the interval ( 1251 ) of 
9— junction existence subdivides into two parts. With a fur- 
ther decrease of £ these parts are transformed into narrow 
bands, which are localized in the vicinity of the — % tran- 
sition point (An + Ap = 0); they take place at kL = k/4 and 
kL = 5n /4. The width of the bands decreases with decrease 
of £. 

Thus, our analysis has shown that for 
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FIG. 3: Analytically derived amplitudes A and B in the CPR of ramp 
S-NF-S structure (d^ = 0.1^/v, dp = 0.65^) and their components 
A/v, Af, Afw versus electrode spacing L at T = 0.77c. Also en- 
hanced interval of <p-state, AL, is marked. 



0.12 <£ 5 0.2 



(26) 



we can expect the formation of 9— junction in a sufficiently 
wide range of distances AL between the electrodes determined 
by (l23l . Now we will take into the account the impact of the 
interface term Ifn{<P)- In the considered approximations, it 
follows from dBT31 )- dBTTb that 



Ifn2(<P) - 



2{/ <?Fexp( 



"-5- J 



I kL 



YBFYBN^I 2 



sin(<p) : 



V2U & . , v „, . (P , 
4^I 2 YbnYbfn^n KVJ V 2 h 



hm{<p) 



2U exp(-f) 



hjBF 



■sin(<p)/sT(sin~), 



(27) 



(28) 



(29) 



where Uq — AW /eyisFNPF- In the range of distances be- 
tween the electrodes 7t/4 < kL < 5n/4 currents Ifni^) an d 
Ifni (<p) are negative. These contributions have the same form 
of CPR as it is for the In{<P) term, and due to negative sign 
suppress the magnitude of supercurrent across the junction 
thus making the inequality ( f23l easier to perform. The re- 
quirement B < imposes additional restriction on the value 
of the suppression parameter Jbfn 



Jbfn > 



Pn^n 
hd N p f 



gg Ybn 

^NjBFNh 1 ! 2 JBF 



(30) 



In derivation of this inequality we have used the fact that in 
the range of distances between the electrodes n/4 < kL < 
5tz/4 depending on kL factor in (1291 is of the order of unity. 
It follows from ( f30b that for a fixed value of Jbfn domain 
of ^-junction existence extends with increase of thickness of 



normal films djv and this domain disappears if becomes 
smaller than the critical value, d^c, 



d NC — 



PnI 



hpFjBFN \£,NjBFNh l l 2 



Ybn 
Ybf 



^ . (31) 



The existence of the critical thickness d^ follows from the 
fact that the amplitude B in In is proportional to d^, while 
in I fn term the parameter B is independent on d^. The sign 
of Ifn\{<P) is positive for n/4 < kL < 3n/4 and negative 
for 2>n/4 < kL < 5n/4 thus providing an advantage for a 
9— junction realization for the lengths which belong to the 
second interval. 

Figure [3] illustrates our analysis. The solid line in Figf3]is 
the modulus of the amplitude of the first harmonic in CPR as 
a function of distance L between S electrodes. It is the re- 
sult of summation of the two contributions following from 
Eqs. ( fTTI i (dash-dotted line) and ( TT2l (dashed line). The 
dash-dot-dotted line in Fig [3] is the amplitude of the second 
harmonic of the CPR following from ( fT2l . The dotted line 
is I F n(L) calculated from flD, ( iBTBI l- ( TbT7T i. All calcula- 
tions have been done for a set of parameters d^ = 0.1 
d F = 0.65^, Jbn = 0.1, Ybf = 1, Ybnf = 10, t, F = O.l^v, 
Pn = Pf, T = 0.77c, H = lOTc- These parameters are close to 
those in real experimental situation. All the amplitudes were 
normalized on factor (lep^ / (WTc))~ l ■ It is evident that there 
is an interval of L, for which the currents in N and F layers 
flow in opposite directions. As a result of the addition of these 
currents the points of — % transitions start to be closer to each 
other. It is seen that in the entire region between these points, 
the inequality (0 is fulfilled. This is exactly the L— interval, 
inside which a <p— junction can be realized. It is also seen that 
contribution of If n part into the full current is small and in 
accordance with our analisys does not play a noticeable role. 

The boundary problem j4l-(|9bl has been solved numeri- 
cally for the same set of junction parameters except c/f . The 
results of calculations for dp — 1 .06^ and dp — 1 .4^ are 
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FIG. 4: Numerically calculated amplitudes A and B in the CPR of 
ramp S-NF-S structure (dfj = O.I^at, dp = 1.06^/v) and their com- 
ponents Aff, Af, Bff, Bf versus electrode spacing L at T = Q.ITq. In 
correspondence with Fig[3]parameters are chosen to form enhanced 
<p-state interval marked by "AL". 



shown in Fig0] and Fig|5] The solid lines in Fig0] are the 
modulus of the amplitudes of the first, A, and the second, 
B, harmonic of CPR as a function of distance L between S 
electrodes. The dashed and dash-dotted lines demonstrate the 
contributions to these amplitudes from the currents flowing 
in N and F films, respectively. All the amplitudes were nor- 
malized on the same factor (2ep^ / (WTc)) . It is seen that 
the main difference between analytical solutions presented in 
Fig|3]and the curves calculated numerically are located in re- 
gion of small L. It is also seen that amplitudes of first and sec- 
ond harmonics of the part of the current flowing in the N film 
slightly decay with L increase. The points of — % transition 
of the first harmonic amplitude of the part of the current flow- 
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FIG. 5: Numerically calculated CPR amplitudes A and B versus elec- 
trode spacing L for S-FN-S structures with dp = 1.06£jv (solid and 
dashed lines respectively) and df = 1.4§y (dash-dotted and dotted 
lines). It is clear that enhanced (p-interval AL\ formed in the first 
case is much larger than pair of ordinary <p-intervals AL2 and AL3 in 
the second one. 



ing in the F layer is slightly shifted to the right, toward larger 
L. It is also seen that the amplitude of the second harmonic, 
Bp, in the interval of interest in the vicinity of L m Q.2^n is 
negligibly small compared to the magnitude of, Bpj. As a re- 
sult, the shape of A(L) curves in Fig J5] and Fig 0] is nearly the 
same, with a little bit larger interval of <p— junction existence 
for the curve calculated numerically. 

Figure [5] demonstrates the same A(L) and B(L) dependen- 
cies as in Fig|4] (solid and dashed lines) together with A{L) 
and B(L) curves calculated for dF = 1.4^ (dash-doted and 
dotted lines). It is clearly seen that for larger dF we get out of 
the interval (l26b and instead of relatively large zone AZ4 may 
have <p— junction in two very narrow intervals AL2 and AL3 lo- 
cated in the vicinity of — % transitions of the first harmonic 
amplitude A . 



V. RAMP TYPE OVERLAP (RTO) JUNCTIONS 

Conditions for the existence of <p— junction d25l ). d26i i can 
be improved by slight modifications of contact geometry, 
namely, by using a combination of ramp and overlap config- 
urations, as it is shown in FigJTJ). Fig|6]demonstrates numer- 
ically calculated spatial distribution of supercurrent in RTO 
^-junction at Josephson phase (p = 7C/2. The current density 
is presented by darkness and the arrows give flows directions. 
The relative smallness of the first harmonics amplitude is pro- 
vided by opposite currents in N and F films. The main fea- 
ture of the ramp-overlap geometry is seen to be specific cur- 
rent distribution in the normal layer leading to another CPR 
shape with dependence on thickness d^. Further, the current 
In should saturate as a function of d^, since normal film re- 
gions located at distances larger than %^ from SN interface are 
practically excluded from the process of supercurrent transfer 
due to exponential decay of proximity-induced superconduct- 
ing correlations 71 . The specific geometry of the RTO struc- 
tures makes theoretical analysis of the processes more com- 
plex than in ramp contact. Nevertheless, it is possible to find 
analytical expressions for supercurrent in these structures and 
to show that the range of parameters providing the existence 
of <p— state is broader than in the ramp type configuration. 

To prove this statement, we consider the RTO structure in 
most practical case of thin N film 

d N « (32) 

and sufficiently large Ybfn providing negligibly small sup- 
pression of superconductivity in N film due to proximity with 
F layer. We will assume additionally that electrode spacing L 
is also small 

L « §v, (33) 

in order to have nonsinusoidal CPR. Under these conditions 
we can at the first step consider the Josephson effect in over- 
lap SN-N-NS structure. Then, at the second step we will use 
the obtained solutions to calculate supercurrent flowing across 
the F part of the RTO structure. The details of calculations are 
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FIG. 6: Current distribution along RTO-type SN-FN-NS structure at 
L = 0.63<^at, dp/ = dp = and T = 0.77c. The intensity of 
gray color shows current density in direction indicated by arrows. 




FIG. 8: The amplitudes of CPR harmonics A, A^, Af. B versus 
electrode spacing L for RTO structure at T <C Tc, Ybm = 0.64 and 
e = 0.123. The mark "AL" shows enhanced <p-state interval. 



summarized in AppendiceslClandlDl They give that the super- 
current 



Is((P)=In((P)+If((P)+Ifn(<P) 



(34) 



consists of three components. Expression for the part of cur- 
rent flowing across N film has the form 



2eI N (q>) 
nTWd N 



r 2 8 2 sin (py/(&YBM + G s ) 



(35) 

jBNdN/^N and jj. = 



w here r = G s / {£IJbm + G s ) , Ybm 
y/£l 2 + r 2 8 2 cos 2 ((p/2), 8 = A/kT c . 

The If((P) term in (f34b is the current through one dimen- 
sional double barrier SFS structure defined by Eq. ( fTTI i, while 
Ifn(<P) is FN-interface term shown in|D] We provide sufficient 
smallness and neglect it in the following estimations. 

As we discussed above, the larger the relative amplitude of 
the second harmonic (or the lower the temperature of a junc- 
tion compare to Tc), the better the conditions for the imple- 
mentation of a (p-contact. In the limit T -C 7c we can trans- 
form from summation to integration over CO in d35l l and calcu- 




-0.2 I 



1 2 3 Jbm 4 

FIG. 7: The amplitudes of the first harmonic Y^ (solid line) and the 
second one Y# (dashed line) normalized on 2W A/ epwYBN versus re- 
duced thickness Ybm- Inset shows the ratio of harmonics [Tj/Xi| 
versus Ybm- 



late numerically the dependence of amplitudes A and B 



epNjBN 



Bn — 



2WA , 
epNjBN 



(36) 



(37) 



on suppression parameter Ybm- The calculated dependencies 
of functions Ta(Ybm) and |Tb|(7bm) are presented in Fig|7j It 
is seen that both Y^ and |Yg| increase with increasing of Ybm 
and saturate at Ybm ~ 1- Inset in Fig|7] shows the ratio of the 
harmonics |Yb/Ya| as a function of Ybm- It achieves maxi- 
mum at Ybm ~ 0.64, thus it determines the optimal values of 
normalized amplitudes of the first Y^ m 0.844 and the sec- 
ond Yg Ri —0.175 harmonics of the current flowing in the N 
layer. It is seen from the inset in Fig|71 that the ratio |Yb/Y a | 
is slowly decreasing function of Ybm- Therefore, the estimates 
given below for Ybm — 0-64 are applicable in a wide parameter 
range 0.5 < Ybm < 10. 

Taking into account these values, we can write down the 
condition of <p-state existence similar to d23l 



e 



<2\Y E 



VfryJf &pf 

Ybn d F p N ' 



(38) 



S?(L) =exp(-fcL)cos^)fL+^ , 

with slightly modified dimensionless parameter e. The wide 
region of <p-state still exists if e is within the interval 

0.123 <£ < 0.298 (39) 

for kL that satisfies the condition ( 1381 . As follows from d38l . 
interval of kL product gains its maximum length 

0.94 < kL < 2.72, (40) 

at e = 0.123. It is seen that these intervals are slightly larger 
than those given by ( |25] | for the ramp type geometry. 
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Fig[8] shows the interval of <p-state existence, AL, in the 
ideal case of T < T c , Jbm = 0.64 and e = 0.123. The cor- 
responding set of parameters dtj = 0.64^, dp = 1.45^, 
Ybn = 1, 7bf = L £f = 0.1^, p w = Pf,H = 10r c was sub- 
stituted in ( fTTI i. d35l >. The solid line is a modulus of the first 
harmonic amplitude, A, its normal, An, and ferromagnetic, 

, parts are presented by dashed and dash-dotted lines re- 
spectively. Finally, the second harmonic amplitude is shown 
as dash-dot-dotted line. It's clear that |A| is relatively small in 
the wide region AL and reaches the value of \2B\ only at local 
maximum. The increased width of AL (see Eqs. (29), (49)) is 
provided by geometric attributes of RTO type structure. 

Let us illustrate the range of nontrivial ground phase (p g ex- 
istence in the structure described in FigJH] The total supercur- 
rent Is is shown on Fig|9]as a function of Josephson phase <p 
and electrode spacing L. It means that each L— section of this 
3D graph is CPR. Solid lines mark the ground state phases at 
each L. In the range of small and large spacing L ground phase 
is located at <p g — 0. However, in the AL-interval CPR be- 
comes significantly nonsinusoidal and demands ground phase 
(p g to split and go to n from both sides; then it— state is real- 
ized at kL = n/2. Clearly, for e > 0.123 the value <p g = n can 
not be reached (see Fig(9^), while in the case of e < 0.123 the 
prolonged 7C-state region is formed (see Fig(9p). 

VI. DISCUSSION 

We have shown that stable <p-state can be realized in S-NF- 
S structures with longitudinally oriented NF-bilayers (though 
<p-state can not be achieved in conventional SNS and SFS 
structures). We have discussed the conditions for realization 
of p-state in ramp-type S-NF-S and RTO-type SN-FN-NS ge- 
ometries. 

Let us discuss most favorable conditions for for experi- 
mental realization of ^-junction. We suggest to use Copper 
as a normal film (<!j# f=s 100 nm and p = 5* 10~ 8 £tm) and 
strongly diluted ferromagnet like FePd or CuNi alloy (£,f sa 
10 nm, H « 10L C ) as the F-layer. We chose Nb (T c « 9K) 
as a superconducting electrode material since it is commonly 
used in superconducting circuits applications. We also pro- 
pose to use sufficiently thick normal layer, above the satura- 
tion threshold, when N-layer thickness have almost no effect. 
After substitution of relevant values into d39l and ( |40l we ar- 
rived at a fairly broad geometrical margins, within which there 
is a possibility for implementation of ^-junctions 

dfii > 50 nm, 
60 nm < d F < 150 nm, (41) 
7 nm < L < 22 nm. 

Finally, the last out-of-plane geometrical scale is set as W = 
140 nm. This value maximizes current and conserves the scale 
of structure in a range of 100 nm. The magnitude of critical 
supercurrent in the <p-state is determined by the second har- 
monic amplitude B 

2WA 

Ic~B N = Y B a 1 mA. (42) 

epNjBN 




FIG. 9: The full current Ig versus Josephson phase <p and electrode 
spacing L for RTO structure at T <C Tc, Ybm = 0.64 and at different F- 
layer thickness parameters a) £ = 0. 137, b) e = 0. 123, c) e = 0. 1 1 1. 
The lines mark the ground states phase (p g . 



The spreads of geometrical scales as well as the magnitude of 
critical current are large enough to be realized experimentally. 

By creating <p-state in a Josephson junction one can fix cer- 
tain value of ground phase (p g . Temperature variation slightly 
shifts the interval of relevant 0-7C transition and permits one 
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to tune the desired ground state phase. Furthermore, sensi- 
tivity of the ground state to an electron distribution function 
permits <p-junctions to be applied as small-scale self-biasing 
one-photon detectors. Moreover, quantum double- well poten- 
tial is formed at the point of ground state splitting providing 
necessary condition for quantum bits and quantum detectors. 
To summarize, Josephson ^-junctions can be realized using 
up-to-date technology and may become important basic ele- 
ment in superconducting electronics. 
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dR F 
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dRpj G F 
dz Gn 



9K(r f -®R n ), <A8a) 

Lrp V tB 



Gf f 03 \ 
-^(Rn-^Rf) 

Gn V CO J 



dU N G F CO 

JbnfSn^— = -pr\UN- —Uf 
dz Gn V co 
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dz G F \ co 
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Appendix A: Ramp type junctions. Limit of small L. 

In the limit of small spacing between S electrodes 

L<min{^,§v} (Al) 

we can neglect nongradient terms in 

d ( 2 d \ d f 2 d \ 
Tx [ Gf ^ Rfn + dz [ Gfn Tz Rf - n = °' (A2) 



Tx{*»Tx Ur * + Tz\r*^z Ur " ]=0 > (A3) 



and introduce four functions 

$> F = R F + iU F , <P n =Rn + iU N , 



(A4) 



where, i, is imaginary unit, R F and Rn are even function of 
coordinate x, while U F and Un are odd in x. Due to the sym- 
metry at x = 



dR 



F.N 



■ o, U FN = 



(A5) 



From ( |A5t and ( |A6at - ( |A7bt it follows that for Jbf and Ybn 
within the interval 



L 



L 



(A10) 



we can neglect Un.f m left hand side of ( IA7al ). ( IA7bb . More- 
over, in this approximation for any point inside the weak link 
region R F ^ 3> U Fi n and the boundary problem (IA2b -( [A9bb for 
functions R F and Rn can be solved resulting in 



R N = Acos(<p/2), R F = — Acos(<p/2) (All) 
CO 



and 



Gn — G F = 



co 



y/co 2 + A 2 cos 2 (cp/ 2) 



(A12) 



Therefore under conditions (lAlOb both Gn and G F are inde- 
pendent on coordinate x,z functions and equations for U f .n 
transform to Laplas equations, which have the solutions 



Asin((j>/2) G s x 
7bn Gn ^ n 

7 -^cosh ?r(2 " +1 L 



U N 

+ £ a n sin ^H+Dl cosh M^ln^^ 



(A13) 



U F 



Asin(ip/2) m Gs x 



Ybf <o Gf sf 
.1 £ ^^(MlScosh!!^. (A14) 



for any coordinate z, and it is convenient to rewrite boundary 
conditions d9al l, d9bl at x = L/2 in the form 

dR G 

Ybn^n~t~~ = 7^ (Acos(<p/2) - R N ) , (A6a) 



Ybf^i 



dx Gn 
dR F Gs ( co 



dx G F \CO 



Acos(<p/2)-fl/r , (A6b) 



Jbn^n^- = ^ (Asin(p/2) - U N ) , 
dx Gn 

. dU F G s ( CO , , . 

^ ar = ^U (<p/ )_ F 



(A7a) 
(A7b) 



They automatically satisfy the boundary conditions at z = 
and z = dN + dp, as well as at x = and x = L/2. To find the 
integration constants a n and b n we have to substitute ( IA1 3b 
and (IA14t into ( |A9at , ( |A9bt and get 



„ _ Asin((p/2)Gs0yjFA/^Ff„ _ t . 7r(2«+l)£4v 
fl "" G^cosh"' 2 "!'^ ' f "- mnn L > 
, _ Asm(<p/2)Gs&7BNF^Ntf . _ . , n{2n+l)d F 



(A15) 



where 



o i k w(2n + l) , , 

P = I YBNFSN p fn + 1 YBFNiFtf + YBNFiNtn, 
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and 



© = 



/ 1 



1 \ 4L (-1)" 



\ JBN^N JBF^F J 7t 2 (2/1 + 1 ) 2 



Substitution of ( IA1 3b and ( IA14I ) into expression for the super- 
current © gives that contributions to the supercurrent across 
the junction proportional to a n and b„ cancel each other and 
Is(<p) equals to the sum 



/s(<p)=/;v(<p)+M<p), 
1 A 2 G N G s sin((p) 



2eI N (<p) = 1 
nTWd N YbnZnPn , 



a>- 



2eI F (<p) = 1 

7iTWd F Ybf^fPf 0)2 



A 2 G F G s sm{(p) 



(A16) 



(A17) 



of the currents, /#(<?>), and, Ip((p), flowing independently 
across F and N parts of the weak link. 



Appendix B: Ramp type junctions. Limit of intermediate L. 

For intermediate values of spacing between the S electrodes 

<^«L«§v. (Bl) 

and suppression parameters at SN and SF interfaces belonging 
to the interval (01 the boundary problem (f4b-d9bb can be also 
solved analytically for sufficiently large suppression param- 
eter Jbfn- Under these restrictions in the first approximation 
we can neglect the suppression of superconductivity in the N 
film due to proximity with the F layer and use expressions 
( IA1 U and ( IA14l i with a„ = as the solution in the N part of 
the weak link. 

To find Rf and Up we have to solve the linear equations 

d 2 d 2 
^-^Rf + ^-^Rf-CIR f = 0, (B2) 



with the boundary conditions 



U F - L~IU F = 0, 



E dR F „ ii A , 
YbfSf^- = G 5 -Acos(<p/2), 
ox Ll 



, dU F „ & t . , ,„s 



at x = L/2, < z < dp and 

dR F 



Ybfn^f 



(B3) 



(B4) 



(B5) 



(B6) 



YBFN%i 



dUp 



= £1 GnUn > 



(B7) 



at z = dp, < x < L/2; (Q = (o/7:T c , £2 = cosign((o) / KT C ) . 
The boundary problem (IB2b - dFJTb must be closed by the con- 
ditions <[Vj and ( IA51 > at free interface of the F film and at the 
line of junction symmetry, respectively. 

Spatial distribution of even in coordinate x part of <P F (x,z) 
can be found in the form of superposition of superconduct- 
ing correlations induced into F film from superconductors and 
from the N part of weak link 



^ _ vnGsAcos((j>/2) 



cosh 



/aG w Acos((j)/2) 



sinh 
cosh 



<£i4- 



(B8) 



&7bfn 



sinh 



Solution for the odd part of <P F (x,z) consists of three terms 



Up 



nV 2 G 6 -Asin((j>/2)^ 

®VBN%NYBFNdF 



, , . xcoshf 1 
VnG s Asin((p/2) V 


^Tf) 


©V sinh I 





-i)"cos 



K-Scoshfic,,^ 



(B9) 



+ 



Vq.i 



Q.G S Asin{<p/2) 



sinh 



^F 



&7bf 



cosh 



where k 2 = L~l + (nn^p/dp) 2 . The first two give the part of 
Up induced from the N film, while the last has the well known 
for SFS junction formed. 

From dB8l and (|B9| | it follows that R*_ m F = R m ,F and 
U*_ mF = U(o, F . Substitution of dB8b and ( IB9b into expression 
for the supercurrent © gives that the Is((p) dependence is 
consists of three terms 



k(<P) = In((P)+If{(P)+Ifn((P)- 



(BIO) 



The first is the supercurrent across the N layer. In consid- 
ered approximation it coincides with the expression given by 
(!Al6l l. The second term in (IB l Oi l is the supercurrent across 
SFS double barrier structure in the limit of small transparen- 
cies of SF interfaces—^ 

2eI F ((p) _ A 2 sin((p) £ G 2 (fin) 



7tTWd F Ybf&Pf 



' or 



£2 sinh (2q F ) 

and the last consists of two terms, I F n{<p) = h(<P) +h{ ( p) 
having different <p— dependence 



2ehM 
nTWd F 



A sin(ip) 



Pf^f 7bf7bfn7bn£n 



(Bl2) 



¥l = 



'a 



2Q. 



sinh(^ L ) sinh(2?i) : 



2el 2 (<p) _ A 2 sin(y) 
nTWd F ~ 7bfnPf(If 



(0= — c 



, o 2 a 2 



^2 



il F n(2q d +smh(2q d )) 
4(/ rf sinh 2 ) 



( l Vp I 

\7bn7bfnZn 



n 



c/ rf cosh(</ i ; 



- 1 



y B F cosh<? L 
2&!;f 



= l Id 



K$ cosh 



( Lk„ \ ' 



(B13) 
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where q^—dpy Q./£f, 1l = Lv Q./2^p. In real experimental 
situation 



The arising boundary problem (lC31 l. iC5\ , iC6\ is also satis- 
fied by independent on x constants leading to 



(B14) 



Under the conditions (IB 14l > some terms of I F n{<?) can be ne- 
glected. Still existing expressions of it parts I F ni (<f>) - Ip^j (<p) 
simplify to 



2eI FN1 ((p) 

TtTWdp 



A 2 sin(<p) %p_ ~ G| Vfl ^ (B15) 
YBFYBFNYBNPFdF %n (o 2 Q. 2 sinhgz.' 



2eI FN2 {(p) ^ A 2 sin((p) y, G N G S % F 
TtTWdp 2y BN Y 2 BFN p F d F w tL„ (0 2 hV 2 <&v ' 



2eI Fm ((p)_ A 2 sin(<p) 



G/vGs 1 



^rWlif YBFNlBFPFdF m ti„ (0 2 Q, CO&hq L ' 



(B16) 



(B17) 



Appendix C: Overlap SN-N-NS junctions. 

To calculate critical current of SN-N-NS junctions we con- 
sider the most practical case of thin N film 



<1n < £a 



(CI) 



and sufficiently large Jbfn providing the absence of suppres- 
sion of superconductivity in N film due to proximity with F 
layer. We will also assume that electrode spacing L is also 
small 



(C2) 



in order to have nonsinusoidal CPR. 

Condition (IC1I ) permits to perform averaging of Usadel 
equations in z— direction in N film, as it was described in de- 
tail in 44 , and reduce the problem to the solution of one dimen- 
sional equations for <$>n — Rjy + iUn. The real part of <Pn is 
the solution of the boundary problem 



ZnYbm 



G N (Q.y BM + Gs) dx 



£IG N dx V N ' dx 



dR N 
'~dx~ 



3Ra 



(p L 

-fl;v = -rAcos-, -<x<°°, 
(C3) 



= 0, 0<x< 



dR* 



= 0, jc = 0, x-H 



(C4) 



(C5) 



where r = G s / (£2y BM + G s ) , Ybm = YbnAn/^n, 8 = A/ nT c . 

From dC4l i. dC3T l it follows that at < x < L/2 functions R N 
are independent on x constants resulting in 



dR 



dx V 2 



(C6) 



R N = rAcos(<p/2), <x < °°. 
Introducing now new functions, 6 

£2 

Un = jUtan0, Gn = — cos0, 



where ji = ^Q? + r 2 8 2 cos 2 ((p/2), we get 



d 2 
dx 2 



A 2 ^-0-sin(e-0) =0, ^<x<°°, 



L 

2 



(C7) 



(C8) 



(C9) 



(CIO) 



(Cll) 



(C12) 



(C13) 



(C14) 



Solution of Eq. (IC9b can be simplified due to existence of the 
first integral 



cos 9 dx 2 ' 2 ' 



de 

9(0) =0, ^=0, x^oo, 



where 



'ISM 



(£27 fiM + G s )V£2 2 + r 2 5 2 ' 



tan0 = 



*sin(<p/2) 



Solution of Eq. ( IClOl l can be easily found 

„ , . 2x , L . L 
9(x) = -9(-),0<x<-. 



2 V dx 



9 +cos(6-<j)) = 1. 



(C15) 



The constant of integration in the right hand side of (IC15I ) 
have been found from the boundary condition (IC lib , which 
demands 6 — >• then x °°. Further integration in ( IC15b for 

L/2 <x<°° gives 



= 0+4 arctan C2 exp 



x-L/2 



(C16) 



where C2 is integration constant, which should be determined 
from the matching conditions atx = L/2. For C2 they give 



(0 +4 arctan (C 2 )) 



2C 2 L 



1+C 2 A' 



(C17) 



Assuming additionally that Ybm is not too small, namely that 
L<S^^ N min ( 1 , ^JYbm) , from dC17b it is easy to get 



,'0 L . 

Co = — tan — sin — 

1 4 4A 2 



(C18) 
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resulting in 



0(jt) = — SHl|, 0<X<-. 



(C19) 



From dC 1 91 > it follows that in weak link region |jc| < L/2 



2x <j> a 

U N = — M sin-, G N = — , 
A A jl 

while under the S electrode, L/2 < x < °° 

Uff = Li tan (0 — 4arctan («)) , 



u = tan -t — 



41 



sin 



(C20) 



(C21) 



Substitution of ( IC7I ). ( IC20I ) into expression (O for the su- 
percurrent in the N channel results in 



2eI N (q>) _ 2 
nTWd N Pn^nJYbm . 



r 2 8 2 sin (p^(Q.y BM + G s ) 

2H^ 2 (vn 2 + r 2 5 2 + Ai) 
(C22) 



with the functions Rn, Gjv, and !/# defined by equations fol- 
lowed from the solution of the boundary problem in the N 
layer described in Appendixlcl 



Rn 



■ rAcos(<p/2), 
U N 



Gn = 



-' Gg x 
Gn 5v ' 



a = 



aAsin'(<p/2)^ 

2^/n 2 +8 2 



(D3) 



(D4) 



r 2 8 2 - 



Substitution of ( ID 1 b - dD4l > into expression © gives that super- 
current across F layer in RTO junction consists of the sum of 
7f(<P) and Ifn(<P), where /f(<p) is the current through one di- 
mensional double barrier SFS structure defined by Eq. ( IB lib , 
while Ifn(<P) = h(<P)+h(<P) has the form 



2eh (<P) _ A 2 sin(y) 

nTWd F p F d F Ybf7bfn^n 

¥1 = 



m=-< 

2Q 



. Q. 2 a>- 



(D5) 



sinh(iy L ) sinh(2^i) ' 



Appendix D: Solution in Ferromagnet Layer of RTO junction. 

Spatial distribution of even and odd in coordinate x parts 
of <Pf(x,z) can be found in the form of superposition of su- 
perconducting correlations induced into F film from supercon- 
ductors and from the N part of weak link. It has the same form 
as in dB8b and i 



Rf = 



, , , coshl v Slr- 
/OGgAcos(g/2) V Sf 



HysF 



sinh 



2 ?F 



coshl V £2-^- 



S-Ybfn 



sinh 



SF 



(Dl) 



2e/ 2 (<p) _ A 2 sin(y) 



1 



Pf^f 7b FN 



\y d F&(2qd+sinh(2q d )) 
- 4q d smh 2 (q d ) 



oo a 2 £l 2 
q d cosh(</i 



, YbfnSn 
I" f 



^2 



y BF coshq L J ' 



1<?rf ^cosh(^j 



(D6) 

Application of conditions (IB 14b allows to neglect some terms 
in/ FA ,(<p) =/FAq(<p)+/FAr2(<p)+4w3(<p) and to simplify re- 
maining terms, leading to the following expressions: 



2eI FNl ((p) _ A 2 sin(<p) t, F " aG 2 s \/H 
%TWd F YbfYbfnPfcIf t, N mfloo (0 2 O. 2 sinhq L ' 



2eI F Nl{<P) _ A 2 ship) ~ raG N G S | F 



(D7) 



(D8) 



coshl A 
VCIGnUn V 



n 3 / 2 g 2 G w (t/ w /x) 

^■Ybfn'If 



sinh 

(-1)" 



SF 



sinh 



Kg cosh K„ 



'Q.G s 8sm(<p/2) 



sinh 



J 



SF 



^rsF 



(D2) 
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